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To develop a sound base on which to design efficient lightguides, 
it is necessary to understand the physics of the propagating modes of 
this transmission medium. In this paper, we examine the properties 
of the scalar modes in an infinite lens- like medium and the vector 
modes in a self-focusing fiber with an infinite homogeneous cladding 
for a hyperbolic secant index distribution. We find that only the 
meridional rays are equalized and that this only occurs when there 
is no material dispersion. 

I. INTRODUCTION 

As is well known, the index distribution 

N(x) = N„sechgx (1) 

has an ideal focusing property for a two-dimensional medium. 1 " 4 The 
group velocities of the various modes are equalized, and if the medium 
is dispersion-free the group velocities are independent of frequency as 
well. It is believed that, for a cylindrical lightguide where the index is 
a function of the radius, such an optimum distribution does not exist. 1 
In particular, a distribution of the sort 

N(r) = No sech gr, (2) 

where r is the radius, will not be expected to have ideal focusing 
properties. In fact, in the presence of material dispersion, the distri- 
bution described by eq. (2) is expected to be similar to a parabolic 
distribution. 

The properties of the distribution described by eq. (2) can most 
easily be inferred from an analysis of the scalar modes in an infinite 
lens-like medium. 1,5 " 9 The purpose of this work is to review these 
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properties and then compare them with an exact numerical analysis of 
the vector modes in a self-focusing fiber with an infinite homogeneous 
cladding. 

II. RELATIONSHIP BETWEEN SCALAR AND VECTOR MODES 

The vector modes 10 " 12 in a lightguide are normally specified as HE m , 9 , 
EH m , 9 , TEo.g, and TMo,,,, where m is the polar index and q the radial 
index. Far from cutoff, relationships between the radial field functions 
for the vector modes </> and the scalar modes \p can be established. 13 
They are 

HE m , 9 modes 

<|>m, 9 — \pm-l,q-l (3) 

EH m , 9 modes 

<t>m,q — ^>m+\,q-l (4) 

te 0i(? modes 

<l>o, Q » ^1,9-1 (5) 

TMo,9 modes 

4>o,q « \Jtl0-l. (6) 

Likewise, relationships can also be established between the effective 
indices N e . They are 





Vector Modes 


Scalar Modes 




HE 


N e (m, q) 


N e (m - 1, q - 1) 


(7) 


EH 


Ne(m, q) 


N e {m + 1, q - 1) 


(8) 


TE 


N e (o, q) 


N e (l, q-D 


(9) 


TM 


N e (o, q) 


Ned, q - 1). 


(10) 



III. SOME PROPERTIES OF SCALAR MODES 

Consider now the following dispersion-free index distribution for an 
infinite lens-like medium 

N 2 (r) = 2V|[1 - (r/L) 2 + 8(r/L) 4 ], (11) 

where 2ttL represents the paraxial focal length of the medium. When 
5 = %, this function approximates eq. (2) closely and, when 8 = o, it 
approximates a parabolic index distribution. For the scalar radial field 
function i^, one finds from a first-order perturbation calculation that 
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the effective indices can be written as 1 

*-*-^ + *5ff- (12) 

where A = 2(2/8 + a + l)/L, (13) 

B = 8[6P 2 + 6jff(a + 1) + (a + l)(o + 2)]/L 2 , (14) 

and A is the wavelength. If we substitute 8 = % and let a = 1, we get 
the result 

This tells us that the difference in effective index between adjacent 
scalar modes with a = 1 is a constant. By means of (9) and (10), we 
predict the identical result for the tightly bound te and tm vector 

modes. 

Also, from first-order perturbation theory, it can be shown that 

N [(2tt/\).No-A/2] 

Ng ' [(2v/\) 2 .N* - (2tt/\)N A + B] 1/2 ' 

where N e is the group index. Again, if we substitute 8 = % and let 
a = 1, we obtain another interesting result, namely, 

N g = N . (17) 

This tells us that all scalar modes with a = 1 have the same group 
index and this index is N„. Equations (5) and (6) predict the same 
result for the tightly bound te and tm vector modes. One would not 
expect this to be the case for general eh or he vector modes, however. 
This latter point will be examined by the numerical methods to be 
described shortly. As pointed out by Kawakami, 1 scalar waves ip a , (l 
when a as o are meridional rays and when /? » o are helical rays. From 
eqs. (3) to (10) and (17), we conclude that the tightly bound te and tm 
vector models are meridional rays. Similarly, tightly bound eh and he 
vector modes with high polar indices m and low radial indices q would 

be helical rays. 

Let us now apply the previous results for scalar waves to a lightguide. 
We consider a guide with a core center index of 1.47428 and at a 
wavelength of 0.82 microns. The parameters N„ and L in eq. (11) are 
chosen so that, at an r of zero, we obtain the core center index and at 
an r of 25 microns we have an index of 1.45330. This allows us to make 
a good comparsion with the vector modes in a self-focusing fiber that 
is considered in the next section. We again assume the medium is 
dispersion-free. 

Figure 1 shows the group index N e of helical and meridional rays as 
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Fig. 1 — The group index N„asa function of effective index N e as calculated from eqs. 
(12) and (16) for meridional and helical rays. The meridional rays are plotted as o while 
the helical rays are plotted as*. In (a), 5 = and the profile is parabolic. In (b), S = Mi 
and there is a separation between helical and meridional rays. In (c), S = % and the 
profile is a hyperbolic secant. We note that the group index of the meridional rays is 
constant. In (d), S = 1 and we note that the helical rays have the same group index. 



a function of effective index for four values of 8, namely, 0, l k, %, and 
1. When 8 = o, we have a parabolic index distribution and we see that 
the helical and meridional rays he on coincident curves. In fact, all 
other rays are coincident as well. We note substantial dispersion of the 
group index with effective index. 

When 8 = V6, the curves for the helical and meridional rays begin to 
separate. The other rays distribute themselves between these two 
curves. One still notes substantial dispersion of the group index with 
effective index. When 8 = %, we very closely approximate a hyperbolic 
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Fig. 1 — continued. 
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secant index profile. The meridional rays now all have the same group 
index or, in other words, are equalized. The helical rays and all other 
rays show varying degrees of dispersion. 

When 5=1, the helical rays are well equalized but now the 
meridional rays and all other rays show dispersion. Quite obviously, 
there is no distribution in index of the form described by eq. (11) that 
would equalize the scalar modes in an infinite lens-like medium. 

III. VECTOR MODES 

Using numerical methods, we can obtain rigorous solutions to the 
vector form of Maxwell's equations for any index distribution. 14-16 The 
method we employ is similar to that of Vigants and Schlesinger 14 and 
uses an optimized fourth-order Runge-Kutta procedure 17 with double 
precision arithmetic. The group indices can be calculated from the 
effective indices by the following formula: 

N g = N e - \dN e /d\. (18) 

We consider the exact hyperbolic secant profile, eq. (2), rather than 
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Fig. 1 — continued. 

the approximate one given by eq. (11) with 8 = %. The wavelength and 
core index are as given previously. However, now we assume a finite 
core with a radius of 25 microns and an infinite, dispersion-free 
cladding with the same index as before. The parameter g in eq. (2) is 

1686 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1980 






chosen so that, when r = 25 microns, the function gives the correct 
value for the cladding index and N is the index at the core center. We 
also consider an exact parabolic-index profile rather than that given 
approximately by eq. (11) with 8 = 0. This profile is tailored to match 
the core-center index and the cladding index. 

From the discussion in the previous section and eq. (15), we would 
expect the difference in effective index to be the same between adjacent 
te or tm vector modes for the hyperbolic secant profile. In Fig. 2, we 
plot this difference for the sequence of tm modes. We see that, for all 
except the last two modes, this spacing is indeed constant. The same 
result is found for the te modes. These latter two modes cannot be 
considered tightly bound because of the finite size of the core of the 
lightguide. Thus, (15) would not apply to them. This change in char- 
acter of these modes clearly demonstrates the strong influence of the 
core-cladding interface. The spacing between adjacent te or tm vector 
modes for the parabolic index profile would not be expected to be the 
same. That this is the case is also shown in Fig. 2. The influence of the 
core-cladding interface is evident here as well. 

From eq. (17) and earlier discussion, we would expect that all the te 
and tm vector modes for the hyperbolic secant profile would have the 
same group index. In Fig. 3, we plot the group index for the te modes 
as a function of effective index. For all practical purposes, the group 
indices are the same except again for the last few modes. The result is 
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Fig. 2 — Plot of AN, versus N e for the te or tm vector modes for a fiber with a 
hyperbolic secant and parabolic index profile, te and tm vector modes correspond to 
meridional rays. This shows that there is no difference between the effective indices of 
adjacent, tightly bound modes. 
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Fig. 3 — Plot of N a versus N e for the TE or TM vector modes for a fiber with a hyperbolic 
secant and a parabolic index profile. For the hyperbolic secant profile, we see that the 
meridional rays have the same group index except for those that are influenced by the 
cladding. 



identical for the tm mode. The average value of group index, if we do 
not include the last modes, is 1.47427. From eq. (17), we would expect 
this to be 1.47428. The agreement is excellent. It should be remembered 
that the TE or tm modes are the meridional rays. Thus the meridional 
rays are seen to be equalized in a lightguide with a hyperbolic secant 
profile, except, of course, the last few, which are being influenced by 
the cladding. Figure 3 also shows the result for the meridional rays 
when the profile is parabolic. Clearly, these rays do not all have the 
same group index. 

From (16), we expect only the te and tm modes to be equalized in 
a fiber with a hyperbolic secant index distribution. In other words, 
skew rays would have different group indices from meridional rays. 
That this is indeed the case can be seen in Fig. 4. The dots correspond 
to the group indices of the vector modes with polar indices 4, 8, 12, 14, 
16, 18, and 20. The solid lines are calculated from eq. (16). The upper 
one is for helical scalar modes, and the lower one for meridional scalar 
modes. It is clear that the vector modes for non-meridional rays are 
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not equalized, and the spread is very similar to that predicted by the 
scalar theory. The strong influence of the cladding on the vector modes 
is again made very clear. 

When material dispersion becomes dominant, such as would be the 
case at 0.82 micron for a germania- doped silica lightguide, we would 
expect great similarity between the group indices of a fiber with a 
parabolic index and a fiber with a hyperbolic secant index distribution. 
Figure 5 shows the group indices for the te vector modes for both 
these profiles with material dispersion effects included exactly in the 
calculation. The material dispersion parameters are taken from the 
work of Fleming 18,19 of 13-mole percent doped silica and pure silica. 
The very strong effect of the materials dispersion is evident. In 
addition, the influence of the cladding on the last two modes is again 
noted. 

IV. CONCLUSIONS 

The scalar modes in an infinite lens-like medium can indeed predict 
some of the properties of the vector modes in a self-focusing fiber with 



1.4743 - 



1.4742 - 



1.4741 - 



1.4740 - 



1.4739 



1.4738 




CORE INDEX 



CLADDING INDEX 



T 



km 
I 

I 
I 



J_ 



1.467 



1.463 

N P 



1.459 



1.455 



Fig. 4 — The dots correspond to the group indices of the vector modes with polar 
indices 4, 8, 12, 14, 16, 18, and 20. The solid lines are scalar modes calculated from eq. 
(16). The upper line is for helical rays and the lower one for meridional rays. 
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Fig. 5 — A comparison of the group indices for a hyperbolic secant profile and a 
parabolic profile in the presence of material dispersion. The wavelength is 0.82 micron, 
and the core center is silica-doped with 13-mole percent germania. 

an infinite homogeneous cladding. Perhaps the most striking difference 
are those modes that are influenced by the cladding. A more subtle 
difference is that the vector modes tend to be more dispersive, as is 
evidenced by Fig. 4. The belief that, for a cylindrical lightguide, a 
distribution of the sort given by eq. (2) (i.e., hyperbolic secant distri- 
bution) only equalizes meridional rays is well substantiated by the 
rigorous vector mode analysis. Finally, we see that in the presence of 
strong material dispersion there is no outstanding advantage of the 
hyperbolic secant-index distribution over the parabolic-index distri- 
bution. 
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